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Preface 


ARCON  Corporation  is  pleased  to  submit  our  finished  Final  Report,  which  describes 
research  that  has  been  carried  out  under  the  provisions  of  contract  F19628-00-C-0072. 
The  period  of  performance  for  this  contract  was  from  July  20,  2000  until  July  24,  2002. 
The  overall  purpose  of  the  work  that  we  have  undertaken  is  to  apply  radiative  transfer 
techniques,  in  new  ways,  to  important  problems  relating  to  background  clutter  mitigation. 
In  particular,  we  simulated  the  radiant  signature  of  gravity  waves  in  the  MWIR  and  used 
our  experience  with  non-LTE  processes  to  investigate  their  likely  contributions  to  the 
observations.  We  also  examined  SABER  limb  sounding  data,  looking  at  the  range  of 
radiance  variability  encountered. 

In  the  process  of  working  with  the  SABER  team,  we  developed  two  new  codes  to  en¬ 
hance  the  capabilities  of  our  non-LTE  model,  Atmospheric  Radiance  Code  (ARC).  These 
are  briefly  described  in  Section  2. 

The  stated  purpose  of  investigating  background  clutter  is  well  served  by  the  work 
reported  below.  It  is  supplemented  by  the  transition  of  some  of  our  modeling  and 
simulation  capabilities  to  the  Air  Force,  in  the  form  of  codes  that  we  have  developed  for 
the  purpose  of  carrying  out  the  work  outlined  above,  and  advice  or  assistance  related  to 
their  use.  Up-to-date  source  code  for  these  applications  has  been  submitted  to  the  Air 
Force  on  a  ZIP  disk.  The  contents  of  the  disk  are  outlined  in  Appendix  A. 


1.  Background  Clutter  in  non-LTE  Regions  of  the  Atmosphere 

Variability  in  atmospheric  infrared  radiance  occurs  on  all  spatial  and  temporal  scales. 
On  the  largest  scales,  seasonal  changes  occur  within  broad  latitudinal  regions.  On  a 
somewhat  smaller  scale,  diurnal  effects  and  planetary  waves  modulate  the  emissions  in 
ways  that  are  at  least  somewhat  predictable.  At  the  same  time  even  smaller  vertical  and 
horizontal  structures  are  often  superposed  on  the  mean  state,  producing  pressure/tem¬ 
perature  perturbations  and  other  localized  effects  that  result  in  substantial  changes  in  the 
emissions  and  the  radiance  levels  observed  by  satellite-borne  remote-sensing  instruments. 
Examples  of  such  structures  are  persistent  temperature  inversion  layers  in  the  meso¬ 
sphere,  gravity  waves,  and  transient  frontal  disturbances.  In  addition,  there  are  bright 
local  phenomena  that  have  little  to  do  with  the  underlying  atmospheric  structure,  the  most 
familiar  example  being  the  aurora. 

The  radiant  signatures  of  these  s mailer-scale  phenomena  offer  a  wealth  of  information 
about  basic  processes  affecting  the  atmosphere,  and  are  often  studied  for  exactly  that 
reason.  At  the  same  time,  they  present  an  obstacle  to  certain  information-gathering  tasks 
that  are  becoming  more  routine  each  year,  because  of  the  unpredictable  background  that 
they  present.  Therefore,  the  need  for  understanding  the  underlying  physics  at  a  very 
detailed  level  is  complemented  by  a  requirement  to  quantify  the  range  and  magnitude  of 
radiance  perturbations  on  a  global  scale. 

We  have  taken  both  the  physical  and  phenomenological  approaches  to  studying  atmos¬ 
pheric  radiance  structure.  In  particular,  we  have  modeled  gravity  waves  in  the  strato¬ 
sphere  and  mesosphere  in  an  attempt  to  quantify  radiance  effects  that  may  be  seen  by 
downlooking  satellite  sensors  operating  in  the  MWIR.  We  also  started  a  project  to  assess 
global  radiance  variability  by  examining  data  from  the  SABER  instrument  on  the  TIMED 
satellite.  These  efforts  are  described  in  the  following  two  sections. 

1.1  Modeling  Gravity  Wave  Signatures 

The  study  of  gravity  waves  was  conducted  in  order  to  simulate  MWIR  radiance  struc¬ 
tures  that  had  been  observed  by  the  MSX  satellite  [Dewan  et  al,  1998\  and  explained  in  a 
preliminary  fashion  by  Picard  et  al  [1998\.  The  data  were  from  MSX  radiometer  band  B, 
which  measured  CO2  4.3  pm  emissions  in  two  slightly  different  but  overlapping  spectral 
ranges.  The  results  of  our  recent  study  have  been  presented  by  Winick  et  al  \2001\  and  by 
Picard  et  al  [2002], 

The  short  summary  of  the  results  can  be  stated  as  follows.  Detailed  line-by-line  calcu¬ 
lations  in  the  MWIR  for  downlooking  (BTH,  or  below  the  horizon)  lines  of  sight  do  re¬ 
produce  radiance  structures  much  like  those  seen  by  MSX.  The  BTH  radiance  response  to 
gravity  waves  in  the  4.3  pm  band  is  largely  determined  by  the  kinetic  temperature  profile, 
insofar  as  most  of  the  emission  seen  by  broadband  sensors  originates  in  the  upper  strato¬ 
sphere  where  LTE  prevails.  However,  a  non-LTE  contribution  in  daytime  from  waves  in 
the  upper  mesosphere  cannot  always  be  ignored,  at  least  for  high  Sun  conditions.  We  also 
showed  that  a  sensor’s  ability  to  discern  wave  structure  is  very  sensitive  to  the  orientation 
of  the  line-of-sight  (LOS)  with  respect  to  the  direction  of  propagation  of  the  wave. 

The  detailed  LOS  calculations  were  performed  using  FASE,  a  line-by-line  (LBL)  code 
incorporating  all  the  line  overlap  and  continuum  effects  needed  for  looking  deep  into  the 
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,  /rACr  .  ,  ,  Contribution  Functions  for  Different  Atmospheric  Profiles 

atmosphere.  (FASE  is  based 
on  the  AFRL  standard,  FAS- 
CODE  [Clough  el  al,  1986], ) 

The  viewing  geometry  was 
chosen  to  mimic  the  BTH 
measurements  with  which  ? 

wave  structure  had  been  seen,  „ 

■a 

namely  those  looking  near  the  « 
horizon  but  still  below  it  at  < 
nadir  angles  of  approximately 
57  degrees.  The  MSX  B1 
spectral  filter  function  [Gar¬ 
diner  et  al,  1994]  was  used  in 
these  calculations. 

Figure  1.  MSX  Bl-band  contribution  functions  for  three  different 

FASE  was  also  used  to  gen-  model  atmospheres, 
erate  contribution  functions 

(CFs).  These  are  quantities  that  record  the  absolute  inband  radiance  contributions  reach¬ 
ing  a  detector  from  within  each  layer  along  a  particular  LOS.  As  such,  they  are  sensitive 
to  the  source  function  at  each  point  along  the  LOS  and  the  optical  path  through  each  la¬ 
yer,  and  their  sum  over  layers  equals  the  total  LOS  radiance.  They  immediately  show  the 
distribution  of  altitudes  within  which  the  detected  emissions  arise.  This  is  how  the  obser¬ 
vations  were  associated  with  upper-stratospheric  thermal  structure  [Picard  et  al,  1998]. 

Figure  1  gives  examples  of  CFs  for  three  model  atmospheres  representing  average 
midlatitude  conditions,  and  polar  summer  and  winter.  The  differences  between  the  curves 
are  largely  due  to  different  stratopause  temperatures;  the  total  radiance  is  much  greater 
for  the  summer  polar  case  than  for  the  others  because  of  this.  The  calculations  were  done 
for  comparison  purposes,  assuming  non-LTE  conditions,  a  nadir  angle  of  57  degrees,  and 
a  solar  zenith  angle  of  60  degrees  (even  for  polar  winter). 

1.1.1  Rapid  Calculations  using  Contribution  Functions 

Contribution  functions  are  also  useful  because  they  enable  us  to  simulate  radiance 
structure  in  a  wave-perturbed  atmosphere  without  repeating  the  tedious  LBL  LOS  calcu¬ 
lations  dozens  or  hundreds  of  times,  as  would  otherwise  be  necessary  for  constructing  an 
image.  The  following  discussion  outlines  the  procedure  for  doing  these  rapid  calcula¬ 
tions. 

It  is  easy  to  show  that,  for  a  LOS  divided  into  N  homogeneous  layers,  the  monochro¬ 
matic  spectral  radiance  is  given  at  the  detector  as  the  following  sum  of  spectral  contribu¬ 
tion  functions  Wvi  [Wintersteiner  and Sharma,  1985], 

W  =  0) 

1=1  i~  1 

where  Rvl  is  the  source  function  in  the  /*  layer,  z w  is  the  optical  path  from  the  near  side  of 
that  layer  to  the  observer,  and  ArV(  is  the  optical  path  within  the  z'th  layer.  The  spectral  con¬ 
tribution  functions  can  be  integrated  over  any  given  spectral  range  Av,  with  a  filter  func¬ 
tion  fv  if  necessary,  to  get  the  total  inband  radiance  as  a  sum  of  layer  contributions,  W, &v,/: 


(2) 


hXs)  =  f  iv(s)fvdv = Z  Ww 

Av  W 

For  LTE  conditions  the  source  function,  Rvi,  is  the  blackbody  function.  As  a  result,  the 
dominant  term  in  W^t  is  nearly  proportional  to  exp(-c2v /7]),  ignoring  the  very  small 

stimulated-emission  factor.  Here,  C2  is  the  second  radiation  constant,  T,  is  the  i  layer 
temperature,  and  v  is  a  representative  energy,  which  we  have  taken  to  be  the  band  center 
at  2349.14  cm'1. 

To  perform  rapid  calculations,  it  is  necessary  to  (1)  generate  CFs,  JFav,/>  numerically 
using  an  unperturbed  or  “background”  model  atmosphere.  FASE  does  this  by  building 
the  total  radiance  from  the  ground  up  at  each  layer  boundary  and  then  subtracting  the 
results  for  consecutive  boundaries,  one  by  one.  [Note  that  it  is  useful  to  have  CFs  for 
more  than  one  nadir  angle  because  this  quantity  varies  along  one  dimension  of  any 
image;  we  did  calculations  for  56.5,  57,  and  57.5  degrees.]  Then  (2)  one  determines  the 
wave-induced  change  in  temperature  and  pressure  with  respect  to  this  background  at 
those  points  where  a  LOS  intersects  each  layer  boundary.  Next  (3)  one  calculates  the 
mftan  source  function  within  each  layer  for  toe  perturbed  conditions,  and  toe  ratio  of  this 
quantity  to  toe  unperturbed  source  fonction.  The  key  step  is  (4)  to  multiply  each  layer’s 
CF  by  this  ratio  and  sum  over  layers  to  produce  toe  perturbed  radiance.  Finally,  (5)  one 
repeats  steps  2-4  for  a  sufficient  number  of  LOS  to  reveal  toe  radiance  structure.  For  a 
single  LOS,  the  perturbed  radiance  so  calculated  is  given  by 

i=i  t=i 

where  f,  is  toe  updated  layer  temperature  reflecting  toe  perturbation  induced  by  toe 
wave. 

Leaving  aside  toe  details  of  step  (2)  for  toe  moment,  some  discussion  of  validation 
tests  we  did  is  called  for.  We  ran  our  tests  with  a  hydrostatically-consistent  midlatitude 
model  atmosphere  that  is  similar  to  toe  U.S.  Standard  Atmosphere,  but,  to  provide  a  more 
realistic  temperature  profile,  lacks  toe  original  discontinuities  in  toe  thermal  gradient.  For 
a  naHir  viewing  angle  of  57  degrees,  toe  simulated  signal  in  toe  MSX  B1  band  from  toe 
frill  LBL  calculation  is  7.769x1 O’7  watt/cm2-sr.  The  first  test  was  to  redo  toe  calculation 
using  toe  contribution  functions,  assuming  no  atmospheric  perturbation.  We  got  a  result 
differing  from  toe  benchmark  by  0.026%  (two  parts  in  toe  fourth  significant  figure).  This 
is  quite  small  compared  to  toe  observed  radiance  structure  amplitude,  which  is  about  3%. 

In  principle,  of  course,  in  toe  absence  of  toe  wave  toe  second  calculation  should  pro¬ 
duce  a  result  identical  to  toe  benchmark.  Possible  sources  of  toe  small  error  we  did  obtain 
include  toe  choices  of  ft  and  v,  and  numerical  roundoff  in  toe  layer-by-layer  subtraction 

that  yields  toe  CFs  in  the  first  place.  We  took  toe  layer  temperature  7]  to  be  toe  average  of 

toe  two  layer-boundary  temperatures  whereas  the  FASE  algorithm  chooses  a  Curtis-God- 
son  weighted  mean  temperature  for  Tt.  With  toe  1-km  layers  specified  in  our  calculation, 
fi-Ti  should  be  very  small.  However,  differences  of  even  a  hundredth  of  a  degree  could 
cause  a  discrepancy  of  ~0.026%  if  they  were  of  the  same  sign  for  most  or  all  of  toe 
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strongly-contributing  layers,  so  it  is  possible  that  revising  the  mean  source  function 
would  reconcile  the  two  calculations.  Use  of  the  mean  band  energy  could  also  skew  the 
result.  A  radiance-weighted  average  over  the  filter  function  might  have  been  a  better 
choice  for  v,  but  it  is  not  obvious  that  any  mean  value  would  lead  to  the  exact  result.  We 
did  not  estimate  the  likely  error  due  to  either  the  choice  of  For  to  roundoff.  Since  the 
numerical  error  is  two  orders  of  magnitude  smaller  than  the  radiance  structure,  we 
proceeded  without  implementing  further  corrections. 

We  also  compared  direct  LBL  calculations  of  the  wave-perturbed  radiance  with  the 
CF-generated  approximation  based  on  the  background  atmosphere,  for  various  phases  of 
the  wave.  These  tests  were  done  assuming  a  wave  with  a  vertical  wavelength  of  10  km 
and  a  horizontal  wavelength  of  32.5  km,  choices  that  were  made  to  maximize  the  observ¬ 
able  contrast  for  viewing  at  the  nadir  angle  of  57  degrees.  The  lower-boundary  forcing 
(see  below)  was  such  that  the  simulated  radiance  structure  amplitude  was  about  0.26x1 0"7 
watt/(cm2-sr),  or  3.3%  of  the  mean,  which  is  comparable  to  the  MSX  observations.  Com¬ 
parisons  between  the  LBL  calculations  and  the  CF  approximation  revealed  differences 
between  0.006xl0'7  and  0.015xl0'7  watt/(cm2-sr),  at  worst  less  than  6%  of  the  amplitude 
or  3%  of  the  peak-to-trough  radiance  contrast.  Thus,  the  CF  algorithm  gives  quite  a 
successful  approximation  to  the  full  LBL  result  despite  the  simple  approach  it  embodies. 
It  is  worth  noting  that  only  the  change  in  atmospheric  temperature  and  not  the  change  in 
density  are  accounted  for  in  these  simulations,  whereas  FASE  accounts  accurately  for  the 
full  effects  of  the  modeled  wave  structure. 

1.1.2  Wave  Model 

Lacking  a  generalized  gravity-wave  propagation  model,  we  adopted  a  simple  algor¬ 
ithm  to  calculate  the  perturbations  in  temperature  and  pressure  that  would  accompany  a 
wave  propagating  through  the  background  atmosphere,  and  then  developed  a  code  called 
GEN_WAVES  to  implement  it.  The  wave  model  is  an  idealization  that  assumes  a  plane 
wave  of  known  vertical  and  horizontal  wavelengths  propagates  in  a  fixed  direction  with 
respect  to  the  LOS.  It  includes  an  arbitrary  forcing  at  the  lower  boundary,  za,  and  assumes 
that  the  amplitude  of  the  temperature  perturbation,  which  is  given  by  A  T(z)  at  any  alti¬ 
tude  z  higher  than  z0,  increases  with  altitude  as  ^N(z0  )  /  N(z)  .  N  is  the  total  number  den¬ 
sity  in  the  background  atmosphere.  The  amplitude  of  the  density  perturbation  is  then 
determined  from  the  low-frequency  approximation  given  by  Walterscheid  et  al  [1987\: 

Atf(z)  =  iV(z)  1-^^j  <♦) 

The  actual  temperature  and  density  perturbations  at  any  point  are  calculated  as  sinusoidal 
functions  (as  shown  below)  of  the  vertical  and  horizontal  wavelengths  and  the  distances 
(x,  z )  from  the  origin.  A  phase  can  be  specified  in  order  to  fix  the  conditions  at  the  origin, 
and  the  forcing  is  adjusted  so  as  to  give  plausible  temperature  perturbations  at  the  strato- 
pause.  It  was  necessary  to  artificially  restrict  the  wave  amplitude  at  high  altitudes  to 
avoid  large  and  unrealistic  temperature/density  swings.  This  was  generally  done  begin¬ 
ning  at  about  80  or  85  km. 

One  important  function  of  the  code  we  wrote  is  to  determine  the  perturbed  atmospher¬ 
ic  state  vector  at  the  intersections  of  the  spherical-shell  layer  boundaries  and  the  LOS,  be- 
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cause  this  is  what  is  needed  to  specify  the  model  atmosphere  for  FASE  or  for  the  approx¬ 
imate  calculation  using  the  CFs.  It  is  an  exercise  in  plane  geometry  to  determine  the  hor¬ 
izontal  distances  from  the  origin  to  these  intersections.  In  the  schematic  drawing  of  Fig¬ 
ure  2,  the  origin  is  chosen  to  be  the  point  of  intersection  of  the  LOS  with  the  Earth’s 
surface.  Measured  from  that  point,  the  horizontal  distance  x„  at  the  rP1  layer  boundary  (at 
altitude  z„ )  is 

Xn  =  Pn{RE+Zn)-  (5) 


In  Eq.(5),  Re  is  the  radius  of  the  Earth  and  P„  is  given  by 


•  ( sin?'  ^ 
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=  arcsin  - - 

~  arcsin 

Ism  Ym) 

V 

Re  sin  Yo 

RE+zn  sin  ym 


(6) 


where  y0  is  the  nadir  viewing  angle  and  ym  is  the  nadir  angle  for  viewing  the  horizon,  e.g. 

'  *  ^  (7) 


ym  =  arcsin 

Hs  being  the  altitude  of  the  satellite. 
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Although  not  indicated  in  Figure  2,  we  also  consider  lines  of  sight  corresponding  to  a 
range  of  nadir  angles,  y.  With  the  origin  still  defined  by  the  intersection  with  the  ground 
of  a  LOS  having  a  nadir  angle  y0,  Eq.(5)  still  applies,  but  the  generalization  of  Eq.(6)  is 


Pn  =  arcsin 
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(8) 


From  this  one  can  calculate  the  wave 
field  at  any  point  (x„,  z„)  on  the  LOS, 
given  the  wavelengths  and  phase.  That 
is,  from  the  altitude-dependent  tem¬ 
perature  and  density  perturbation 
amplitudes  A7\z„)and  AN(zn),  and 
the  assumed  functional  form,  one  can 
easily  derive  the  changes  in  tem¬ 
perature  and  density  and  hence  the 
complete  atmospheric  state  at  the 
relevant  locations.  It  is  true  that  nei¬ 
ther  the  geometry  nor  the  wave  model 
are  general  enough  to  simulate  all 
cases  of  interest,  but  the  purpose  of  in¬ 
vestigating  the  origins  of  observable 
MWIR  radiance  structure  is  served 
nonetheless. 


Figure  2.  Schematic  showing  simplified  geometry  for  BTH  viewing  at  a  nadir  angle  of  y0  from  a  satellite  S. 
Horizontal  and  vertical  coordinates,  with  respect  to  the  origin  at  O,  of  the  intersection  of  the  LOS  with  the 
boundaries  of  the  nth  atmospheric  layer  are  indicated,  as  is  the  angle  8  between  the  LOS  and  the  local  horizon¬ 
tal  at  the  top  of  the  layer. 
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1.1.3  Orientation  of  the  LOS 

Due  the  confluence  of  weak  and  strong  lines,  contribution  functions  calculated  for 
wide  spectral  intervals  display  a  broad  altitude  distribution.  CFs  shown  in  Figure  1  cover 
more  than  20  km  at  half-maximum,  which  is  considerably  greater  than  typical  gravity- 
wave  vertical  wavelengths.  As  a  result,  the  likelihood  of  observing  radiance  structure  due 
to  waves  depends  on  the  viewing  geometry.  When  the  view  intersects  many  wave  peaks 
and  troughs  in  the  vicinity  of  the  CF’s  maximum,  cancellation  occurs.  Then  the  radiance 
is  not  much  different  from  that  of  the  unperturbed  atmosphere  and  the  contrast  between 
lines  of  sight  with  slightly  different  viewing  directions  (due,  for  example,  to  slightly 
different  nadir  angles)  is  minimized.  On  the  other  hand,  LOS  with  just  the  right 
orientation  with  respect  to  the  wave-propagation  direction  may  intersect  only  peaks  or 
only  troughs  over  much  of  the  total  range,  hi  effect,  each  LOS  sees  only  a  single  phase  of 
the  upwelling  wave  at  all  altitudes,  but  the  phase  that  is  seen  is  different  for  slightly  dif¬ 
ferent  viewing  angles.  This  is  the  condition  that  maximizes  the  contrast  between  different 
LOS  and  reveals  the  underlying  property  of  the  atmosphere,  namely  the  temperature 
excursions  that  signify  the  passage  of  a  wave. 

For  a  monochromatic  wave  of  the  sort  we  are  considering  in  our  simple  model,  the 
propagation  direction  is  determined  by  the  ratio  of  the  horizontal  and  vertical  wave¬ 
lengths,  Xx  and  4.  Moreover,  the  spatial  dependence  of  the  perturbation  is  given  by 


f 


AT  ~  sin 


2ZX  2  7TZ 


(9) 


\  "x  J 

where  <f>  is  an  arbitrary  phase  and  the  choice  of  sign  in  the  argument  determines  whether 
the  wave  is  propagating  toward  or  away  from  the  observer.  For  the  phase  of  the  wave  to 
remain  constant  over  a  range  (Ax,  Az)  on  the  locus  of  points  (x,  z)  comprising  the  LOS, 
one  can  only  choose  waves  propagating  toward  the  observer,  e.g.  using  the  negative  sign 
in  Eq.(9).  Then  the  condition  to  be  satisfied  is 


—  «  —  =  cotan(£),  (10) 

4  Az 

where  £is  the  angle  between  the  LOS  and  the  local  horizontal  indicated  in  Figure  2.  For 
the  MSX  geometry,  5  is  about  18.2  degrees  at  40  km  and  the  wavelength  ratio  is  about 
3.2.  This  relation  is  exact  only  in  the  limit  of  small  Az,  since  5  varies  slowly  with  z  due 

to  the  Earth’s  curvature.  But  one  can  also  show,  by  the  tedious  exercise  of  expressing  x  in 
terms  of  z  using  Eqs.(5)  and  (6)  and  equating  the  phases  at  different  altitudes,  that  the 
expression  is  in  error  by  a  very  small  amount  for  values  of  Az  close  to  the  width  of  the 
CFs. 


A  simple  comparison  is  useful  for  illustrating  the  effect  of  the  viewing  direction.  Fig¬ 
ure  3  plots  the  temperature  T  and  the  perturbation  AT  on  the  LOS  near  the  stratopause, 
considering  an  optimally-propagating  wave  as  specified  by  Eq.(10)  for  the  MSX  geome¬ 
try  and  also  a  wave  with  the  same  properties  but  propagating  in  the  opposite  direction.  In 
one  case  the  temperature  on  the  LOS  is  elevated  with  respect  to  the  background  tempera¬ 
ture  at  all  altitudes  in  the  range  shown.  For  slightly  different  nadir  angles  it  would  be  less 
elevated,  or  diminished,  at  all  altitudes.  For  the  case  with  the  wave  propagating  away 
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Midlatitude  Model  with  Wave  Perturbation, 
Profiles  Corresponding  to  LOS 


Perturbation 


Kinetic  Temperature  (K) 


AT  (K) 


Figure  3  (a)  Temperature  profiles  along  the  line  of  sight,  for  the  midlatitude  background  model  and  two 
simulated  waves  propagating  toward  and  away  from  the  observer.  The  calculations  were  done  for  a  nadir 
viewing  angle  of  57  degrees,  assuming  vertical  and  horizontal  wavelengths  of  10  and  32  km  respectively  so 
that  the  phase  fronts  of  one  of  the  waves  would  be  nearly  aligned  with  the  LOS.  (b)  The  perturbation  with 
respect  to  the  background,  to  illustrate  that  the  temperature  is  uniformly  enhanced  on  the  LOS  for  the 
aligned  case  and  not  for  the  nonaligned  case.  See  the  discussion  in  the  text. 


from  the  observer,  however,  the  view  encounters  many  peaks  and  valleys.  Here  the  radi¬ 
ance  is  hardly  different  from  that  of  the  background  atmosphere,  and  slightly  different 
viewing  angles  merely  shift  the  altitudes  on  the  LOS  at  which  the  peaks  appear  without 
appreciably  changing  the  sum  of  the  contributions  from  different  segments  of  the  LOS. 


It  follows  that  the  viewing 
geometry  is  a  critical  consid¬ 
eration  for  detecting  strato¬ 
spheric  waves  through  their 
radiant  signatures.  Winick  et  al 
[2001]  calculated  the  radiance 
as  a  function  of  nadir  angle  for 
the  nearly-aligned  case  depic¬ 
ted  in  Figure  3  and  for  other 
cases  with  different  wave¬ 
length  ratios.  Figure  4,  from 
their  presentation,  shows  the 
radiance  contrast  for  the  near¬ 
ly-aligned  case,  Ax  -  32.5  km, 
A,  =  10  km  on  the  left.  The 


Simulated  Radiance  Across  1  FOV  X  «10km,  Waves  Approaching  Observer 


X  -32.5  X  -25.0  X  -47.0  X  -20.0  background 


Figure  4.  Simulated  radiance  for  a  1 -degree  range  of  nadir  viewing  angles  centered  on  57  degrees,  assuming 
different  horizontal  wavelengths  and  therefore  different  degrees  of  alignment  of  the  phase  fronts  with  the  LOS. 
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cases  with  ^  equal  to  25  and  47  km  (everything  else  remaining  as  before)  correspond  to 
waves  propagating  only  about  5  degrees  from  the  “optimal  direction,  but  even  so  the 
contrast  is  much  reduced.  For  Xx  =  20  km  the  discrepancy  is  about  10  degrees,  and  one 

can  see  that  that  radiance  pattern  is  almost  the  same  as  that  of  the  background  atmos¬ 
phere,  given  on  the  right.  [The  change  in  the  viewing  angle  alone  causes  the  variation 
within  the  rightmost  profile  in  Figure  4.  There  is  a  slightly  greater  airmass  on  the  LOS  for 
larger  nadir  angles,  resulting  in  greater  opacity  and  therefore  CFs  that  peak  at  slightly 
higher  (and  warmer)  altitudes  in  the  stratosphere  ] 

One  can  conclude  that  the  viewing  geometry  imposes  strict  conditions  on  what  can  be 
detected  and  what  can’t.  Therefore,  even  if  upward-propagating  waves  are  ubiquitous  in 
the  atmosphere,  the  probability  of  observing  them  is  much  reduced  by  geometrical/obser- 
vational  constraints. 


1.1.4.  Non-LTE  Effects 

The  discussion  of  contribution  functions,  and  in  particular  their  use  for  rapidly  deter¬ 
mining  LOS  radiance  in  a  wave-perturbed  atmosphere,  omitted  any  mention  of  non-LTE 
effects.  The  relevant  C02  4.3  pm  states  are  in  fact  in  LTE  throughout  the  stratosphere, 
where  the  CF  peaks  occur,  but  in  the  mesosphere  their  vibrational  temperatures  depart 
from  equilibrium  because  of  radiative  transfer  from  the  warm  stratopause  and  because  of 
the  absorption  of  solar  flux.  As  a  result,  the  total  observed  radiance  for  BTH  views  does 
include  some  contribution  from  non-LTE  regions,  although  at  night  it  is  very  small. 


The  non-LTE  influence  can  be  seen  in  the  CF  plots  of  Figure  1,  where  the  secondary 
peak  near  80  km  is  caused  by  greatly-elevated  daytime  vibrational  populations.  (In  fact, 
for  the  polar  winter  model  the  non-LTE  contribution  appears  dominant,  which  would  be 
the  actual  situation  if  daytime  conditions  could  actually  prevail.)  It  is  further  emphasized 
in  Figure  5,  where  CFs  calculated  for  both  LTE  and  non-LTE  conditions  are  compared 
for  the  midlatitude  model  we  use.  The  non-LTE  calculations  assume  that  the  minor- 


isotope  C02  vibrational  tempera¬ 
tures  are  the  same  as  the  major-iso- 
tope  temperature,  a  limitation  of 
FASE  that  introduces  a  small  error. 
But  the  fact  that  there  is  a  signifi¬ 
cant  non-LTE  component  to  the 
total  is  apparent. 

The  implications  that  non-LTE 
effects  have  for  observations  of 
wave  structure  are  not  thoroughly 
understood  at  this  time.  It  is  impor¬ 
tant  to  note,  however,  that  the  kine¬ 
tic  and  vibrational  temperatures  are 
largely  decoupled  above  about  60 
km,  because  the  dominant  excita- 


LTE  vs  NLTE  Contribution  Function:  U.S.  Standard 
MSXB1,  Nadir  Angle  =  57° 


Contribution  Function  (watts-cm'2-sr'/km) 


tion  mechanisms  for  the  principal 
emitting  states  are  not  thermal  as 
they  are  at  lower  altitudes.  As  a 


Figure  5.  LTE  and  non-LTE  contribution  functions  for  the 
midlatitude  model  atmosphere,  from  FASE. 
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Wave-Induced  Perturbations  of  Kinetic  and  Vibrational  Temperatures 


Figure  6.  Wave-induced  perturbations  of  (a)  kinetic  temperature;  (b)  vibrational  temperature,  TV,  of  die 
C02(00011)  state  for  four  isotopes;  and  (c)  Tv  for  the  Group  1  states  of  the  626  isotope.  Most  of  the  non- 
LTE  contribution  to  the  BTH  radiance  arises  in  these  solar-pumped  states  near  80  km.  The  kinetic 
temperature  perturbation  is  included  in  both  (b)  and  (c)  for  comparison  purposes. 


result,  the  wave-induced  excursions  in  temperature,  which  are  much  larger  in  the 
mesosphere  than  below  the  stratopause  (c.f.  Figure  3),  have  relatively  little  effect  upon 
the  emission  rate.  That  is,  even  though  the  total  modeled  radiance  is  greater  because  of 
non-LTE  emissions,  the  non-LTE  component  is  relatively  insensitive  to  temperature 
changes,  and  the  radiance  contrast  between  different  LOS  is  determined  primarily  by  the 
perturbation  of  the  stratosphere  rather  than  that  of  the  mesosphere. 

The  effect  of  the  relatively  weak  coupling  between  the  kinetic  temperature  and  the 
populations  of  the  4.3  pm  states  in  the  mesosphere  is  illustrated  in  Figure  6.  We  calcula¬ 
ted  vibrational  temperatures  of  the  most  important  emitting  states  in  the  presence  of  a 
wave,  assuming  daytime  conditions  with  the  Sun  at  a  zenith  angle  of  60  degrees,  and 
compared  them  with  temperatures  determined  for  the  background  atmosphere.  Although 
the  non-LTE  model  calculation  is  approximate,  the  wave-induced  changes  in  vibrational 
temperature  are  much  smaller  than  the  changes  in  the  kinetic  temperature  above  60  km, 
and  indeed  smaller  than  the  changes  in  vibrational  temperature  at  40  km.  Moreover,  for 
various  reasons,  these  temperature  variations  can  have  different  phases  with  respect  to 
each  other  in  the  mesosphere,  so  cancellation  also  occurs.  Both  of  these  circumstances 
underlie  the  conclusion  that,  despite  the  rapid  growth  of  the  wave  amplitude  in  the  meso¬ 
sphere,  its  effect  on  the  overall  emission  rate  as  integrated  over  the  non-LTE  segment  of  a 
BTH  viewing  path  is  small.  More  to  the  point,  it  is  very  small  compared  to  effects  from 
lower  down  where  the  densities  are  much  greater  and  the  CF  has  its  maximum  value. 

There  are  several  reasons  that  non-LTE  results  are  approximate.  One  is  that  the  hori¬ 
zontal  spatial  variations  introduced  by  the  wave  cannot  be  accommodated  in  the  radia¬ 
tive-transport  calculation  involved  in  determining  the  vibrational  temperatures.  For  the 
non-LTE  region  the  mean  free  path  of  photons  in  some  4.3  pm  bands  may  be  comparable 
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to  the  horizontal  wavelength,  so  in  the  actual  atmosphere  there  would  be  some  averaging 
over  regions  with  different  emission  rates.  Another  is  that  it  is  very  difficult  to  calculate 
the  solar  flux  absorption  rates  while  taking  into  account  the  horizontal  spatial  variability. 
The  calculations  shown  above  assumed  only  vertical  variations.  Perhaps  even  more 
important,  the  non-LTE  signal  depends  sensitively  upon  the  actual  growth  of  the  wave, 
which  is  something  that  can  be  modeled  only  approximately  in  this  region.  Last  of  all, 
FASE  does  not  correctly  treat  the  minor  isotopes  in  the  non-LTE  LOS  calculation,  since 
it  assumes  that  their  vibrational  temperatures  are  the  same  as  that  of  the  major  isotope. 

For  all  of  these  reasons,  it  is  very  difficult  to  accurately  simulate  the  non-LTE  compo¬ 
nent  of  the  BTH  radiance  in  the  presence  of  waves.  We  do  know,  at  least  for  daytime 
conditions  with  a  Sun  angle  of  75  degrees  or  less,  that  the  overall  non-LTE  contribution 
is  a  significant  portion  of  the  total  (as  suggested  by  Figure  5)  and  therefore  must  be 
included  in  some  fashion.  Relying  on  the  conclusion,  stated  above,  about  the  diminished 
wave  influence  in  the  mesosphere,  we  explored  the  possibility  of  combining  the  non-LTE 
contribution  from  the  background  atmosphere  with  the  LTE  contribution  in  the  wave-per¬ 
turbed  atmosphere  to  arrive  at  a  proper  simulated  BTH  radiance.  However,  that  work  is 
not  complete,  and  a  complete  characterization  of  the  non-LTE  effects  remains  a  matter 
for  future  consideration. 

1.2  SABER 

The  phenomenological  approach  that  we  took  toward  studying  atmospheric  radiance 
structure  involved  the  SABER  instrument  on  the  TIMED  satellite.  TIMED  (Thermo¬ 
sphere  Ionosphere  Mesosphere  Energetics  and  Dynamics)  was  launched  in  December 
2001.  SABER  (Sounding  of  the  Atmosphere  using  Broadband  Emission  Radiometry)  is  a 
10-channel  limb-sounding  instrument  whose  overall  purpose  is  to  retrieve  the  temper¬ 
ature/pressure  structure  of  the  atmosphere,  and  to  measure  a  number  of  emissions  that  are 
critical  for  determining  the  energy  balance  of  the  mesosphere/lower  thermosphere  (MLT) 
region.  SABER  provides  nearly  global  coverage  (either  from  ~52°S  to  ~82°N  or  from 
~82°S  to  ~52°N)  on  a  continual  basis.  It  is  developing  a  voluminous  database  of  radiance 
profiles,  which  ultimately  will  be  associated  with  an  even  more  voluminous  set  of 
retrieved  or  derived  quantities.  The  breadth  and  scope  of  the  coverage  that  SABER  is 
providing  is  unprecedented,  for  measurements  of  the  MLT  region. 

One  of  the  many  things  that  will  eventually  be  derived  from  the  SABER  data  is  a 
global  picture  of  the  radiance  and  radiance  variability  in  each  of  its  ten  channels.  The 
variability  will  be  assessable  on  many  of  the  different  scales  we  outlined  at  the  beginning 
of  this  section.  When  it  is  suitably  organized  and  reduced  into  products  that  can  be 
scanned  in  a  systematic  fashion,  it  will  answer  many  outstanding  questions  about  the 
frequency  and  magnitude  of  radiance  perturbations,  and  will  put  limits  upon  what  can  be 
expected  to  be  seen  in  the  future. 

As  soon  as  Level  IB  data  (radiance  data  in  physical  units)  were  available,  which  was 
only  very  recently,  we  began  a  process  of  organizing  the  data  set  and  looking  at  the  radi¬ 
ance  variability.  We  started  with  data  from  a  single  day,  3  March  2002.  The  set  came 
from  all  14  orbits  that  day,  containing  approximately  1360  events.  An  event  is  a  single 
sweep  of  the  sensors  through  a  range  of  tangent  heights  from  hard  Earth  to  ~300-400  km. 
We  looked  at  the  CO2  and  OH  channels.  Channels  1-3,  which  have  different  passbands, 


sense  the  15  pm  emissions  that  are  used  for  (among  other  things)  retrieval  of  kinetic  tem¬ 
perature.  Channel  7  measures  4.3  pm  radiance,  mostly  from  CQ2.  Channels  8  and  9  detect 
Av  =  2  transitions  from  OH  at  2.0  and  1.6  pm,  respectively.  For  channel  8  these  originate 
in  vibrational  states  8  and  9  (e.g.,  the  9-7  and  8-6  Meinel  bands  are  seen)  with  a  small 
contribution  from  v'  =  7  as  well.  For  channel  9  the  bands  are  5-3  and  4-2,  with  some  6-4 
mixed  in. 

We  have  not  looked  at  data  from  channels  4-6  (emissions  from  O3,  H2O,  and  NO 
respectively)  or  at  channel  10,  which  sees  emissions  from  02(1Ag)  at  1.27  pm.  In  the 
ozone  and  water  channels  there  are  known  issues  concerning  off-axis  rejection  that  have 
not  been  resolved,  and  in  any  case  other  researchers  are  looking  at  these  data  in  detail. 
The  NO  data,  and  CO2  4.3  pm  data,  are  of  particularly  high  quality.  Usable  data  extend  to 
much  higher  tangent  heights  than  for  any  previously-accumulated  data  set,  up  to  250  km 
in  the  case  of  NO  and  200  km  for  CO2  in  the  daytime  or  when  viewing  active  aurora. 

The  codes  we  used  to  perform  the  analysis  that  we  present  below,  and  that  are  capable 
of  other  tasks  as  well,  are  included  among  the  software  that  is  being  submitted  to  the  Air 
Force  along  with  this  report.  Appendix  A  gives  a  short  discussion  of  the  tasks  that  these 
codes  perform. 

On  the  next  two  pages  we  show  examples  of  typical  SABER  limb  radiance  scans  for 
several  of  the  channels  we  have  been  evaluating.  Figures  7  and  8  show  nighttime  and 
daytime  events  from  northern  midlatitude  regions.  Figures  9  and  10  are  for  equatorial 
events.  Parameters  that  we  use  to  characterize  each  event— the  SZA,  latitude,  etc— vary 
over  the  extent  of  each  limb  path,  and  are  also  different  for  paths  with  different  tangent 
heights.  The  values  listed  on  the  figures  refer  to  the  tangent  point  of  the  80-km  limb  path 
in  each  case. 

Figure  7  displays  nighttime  data  from  the  four  CO2  channels  and  one  OH  channel  for  a 
midlatitude  event.  The  15  pm  data  on  the  right  mostly  reflect  the  thermal  structure  of  the 
atmosphere.  The  protrusion  that  appears  near  110  km  is  a  common  feature  of  SABER 
scans,  and  indeed  similar  features  have  been  observed  in  past  experiments.  The  signal  to 
noise  ratio  exceeds  10  for  tangent  heights  up  to  about  115  km,  typical  of  SABER  data. 
The  4.3  pm  radiance  on  the  left  is  affected  by  the  kinetic  temperature  profile  and  many 
different  processes,  and  includes  small  iriband  contributions  from  OH  (9-8)  and  (8-7) 
emissions.  Here  the  S/N  exceeds  10  at  120  km.  The  OH  profile  is  also  typical  of  night¬ 
time  events. 

Figure  8  shows  a  daytime  midlatitude  event.  The  channel-7  radiance  is,  as  expected, 
one  to  two  orders  of  magnitude  higher  than  its  nighttime  value  at  all  tangent  heights  in 
the  mesosphere  and  thermosphere.  In  this  case  the  useful  data  extend  nearly  up  to  200 
km.  There  probably  is  a  contribution  from  NO+  emission  above  120  km. 

Figure  9  is  similar  to  Figure  7  but  the  data  were  acquired  near  the  equator.  The  15  pm 
radiances  in  channels  1-3  all  have  a  very  pronounced  layer  structure  at  80  km,  a  feature 
that  was  seen  in  many  low-latitude  events  on  that  day.  The  4.3  pm  and  OH  radiances  also 
have  maxima  at  this  tangent  height;  in  feet,  the  channel-7  radiance  is  a  factor  of  five 
greater  at  80  km  than  it  is  in  the  event  shown  in  Figure  7.  The  110-km  feature  in  the  15 
pm  channels  is  less  pronounced  here  than  in  the  previous  figures,  but  is  still  typical. 
Some  events,  however,  have  virtually  no  change  of  slope  in  this  region. 
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Figure  7.  Typical  midlatitude  nighttime  event,  showing  SABER  radiances  in  CO2  and  OH  channels 
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Figure  8.  Typical  daytime  event,  showing  SABER  radiances  in  the  four  C02  channels. 
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Figure  9.  As  Figure  7,  but  for  an  equatorial  event 
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Figure  10.  As  Figure  8,  but  for  an  equatorial  event 
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Figure  10  shows  a  daytime  equatorial  event.  One  can  see  that  the  layer  appearing  in 
the  upper  mesosphere  in  the  15  |xm  channels  in  Figure  8  is  also  seen  here,  but  at  a  some¬ 
what  higher  tangent  height  (~88  km  instead  of  ~80  km).  Although  not  so  prominent  in  all 
scans  as  in  these  two,  this  structure  appears  in  most  of  the  equatorial  events  that  we 
examined  among  the  March  3  data.  The  present  interpretation  is  that  it  must  be  due  to 
temperature  inversion.  In  feet,  preliminary  LTE  retrievals  of  temperature  reveal  this;  for 
the  event  in  Figure  10  an  increase  of  ~60  K  in  Tt  is  seen  near  88  km.  However,  it  must  be 
noted  that  since  non-LTE  effects  were  not  accounted  for,  all  quantitative  results  like  this 
are  estimates  that  are  likely  to  be  changed  considerably  in  the  future.  Moreover,  day/night 
differences  are  suggestive  of  tidal  effects,  as  discussed  below. 

We  note  that  no  correction  for  off-axis  non-rejected  earth  radiance  has  been  imple¬ 
mented  in  the  version  of  the  radiance  data  shown  here.  For  some  channels  and  some 
ranges  of  tangent  heights,  there  are  significant  contributions  from  scattered  light  origina¬ 
ting  at  low  altitudes.  The  problem  is  most  acute  when  the  radiance  profile  drops  sharply 
with  increasing  tangent  height.  The  magnitude  of  corrections  that  will  be  made  is  not 
known  at  this  time,  but  for  future  data  releases  reductions  in  limb  radiance  of  as  much  as 
50%  are  thought  to  he  possible  for  the  most  strongly-affected  channels.  Since  the  correc¬ 
tions  are  likely  to  be  similar  for  all  similar  events,  the  radiance  variability  that  we  observe 
will  undoubtedly  remain  even  after  they  are  implemented.  However,  the  absolute  radi¬ 
ance  levels  may  finally  be  quite  different  for  some  channels  at  some  tangent  heights. 


SABER  C02  Limb  Radiance 


Figure  1 1.  A  nighttime  high-latitude  event,  showing  evidence  of  auroral  excitation  of  the  4.3  pm  states. 
A  nighttime  simulation  assuming  quiescent  conditions  is  shown  to  illustrate  the  auroral  enhancement. 
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Figure  1 1  shows  a  high-latitude  event  with  a  very  large  enhancement  of  radiance  in 
channel  7,  surely  caused  by  an  aurora.  March  3  was  a  quiet  day  overall,  but  several  high- 
latitude  events  show  elevated  channel-7  radiance,  this  being  the  most  notable.  The  figure 
includes  an  ARC  simulation  of  a  quiet  nighttime  event  to  show  the  size  of  the  effect.  (The 
non-LTE  calculation  for  that  day  gave  excellent  agreement  with  the  15  pm  radiance,  so 
that  result  is  reasonably  reliable.) 

Modeling  the  4.3  gm  signature  of  an  aurora  is  a  very  difficult  task  because  of  uncer¬ 
tainty  about  the  level  of  dosing,  and  the  location  and  extent  of  the  dosing  region  along  the 
LOS.  hi  feet,  a  strong  aurora  on  the  fer  side  of  the  tangent  point  produces  only  a  very 
small  change  in  radiance,  due  to  atmospheric  opacity,  whereas  the  same  dosing  closer  to 
the  observing  instrument  could  cause  an  enhancement  of  an  order  of  magnitude  as  shown 
here,  or  possibly  even  more.  Moreover,  it  is  not  clear  how  much  of  increases  in  radiance 
that  are  observed  are  due  to  directly-imaged  aurora  and  how  much  may  be  due  to  ancil¬ 
lary  effects  from  past  activity,  such  as  atmospheric  heave.  Nevertheless,  it  is  useful  to 
keep  track  of  the  extent  by  which  these  radiance  levels  vary. 

Figure  1 1  also  shows  an  enhancement  in  channel  8  radiance  above  about  100  km.  This 
is  seen  in  other  high-latitude  events  and  correlates  with  elevated  4.3  gm.  It  is  unlikely  that 
OH  is  making  much  of  a  contribution  from  these  altitudes,  but  the  actual  source  of  the 
increase  in  2  gm  radiance  has  not  been  identified. 

Just  from  looking  at  individual  limb-scanning  events,  it  is  clear  that  a  lot  of  variability 
exists  even  among  events  that  one  might  expect  to  be  alike— for  example,  nighttime  events 
at  similar  latitudes,  or  consecutive  events  separated  by  relatively  small  distances.  To 
assess  this  variability,  we  tried  to  devise  ways  of  organizing  and  viewing  data  from  many 
events  at  once.  Figures  12-18  show  one  way  of  doing  this.  In  these  figures,  we  display 
averages,  standard  deviations,  and  minimum  and  maximum  radiances  calculated  from  the 
channel  1  (15  gm)  data  from  all  orbits  of  March  3,  2002,  and  the  same  for  channel  7  (4.3 
gm)  nighttime  data.  Nighttime  events  were  defined  to  be  those  for  which  the  SZA  at  the 
80-km  tangent  point  was  105  degrees  or  greater.  The  figures  break  out  the  events 
according  to  latitude  bins  that  together  cover  the  full  range  of  territory  viewed  by  SABER 
on  that  day. 

Note  that  averaging  channel-7  daytime  data  would  not  be  meaningful  because  of  the 
large  differences  in  the  excitation  of  4.3  gm  states  that  is  expected  for  the  range  of  solar 
zenith  angles  encountered.  The  15  gm  states,  on  the  other  hand,  do  not  respond  directly  to 
solar  influences,  so  to  first  order  one  expects  fall-day  averages  of  channel  1  data  to  be 
useful  quantities.  However,  see  the  discussion  below. 

The  first  few  figures  in  this  compendium  show  data  from  low  and  mid  latitudes.  The 
northern  midlatitude  radiance  (Figure  12)  was  apparently  consistent  over  the  globe,  as  the 
standard  deviations  and  even  the  min-max  envelopes  suggest  rather  small  excursions 
about  the  mean.  This  was  less  true  at  lower  latitudes  on  that  particular  day  (Figures  13- 
15).  In  that  case,  there  were  much  larger  deviations  in  the  upper  mesosphere  in  both 
channels  between  30  degrees  north  and  30  degrees  south,  particularly  in  the  bin  spanning 
the  equator,  and  the  15  gm  radiance  also  has  a  great  deal  more  variability  near  110  km.  In 
the  lower  mesosphere,  on  the  other  hand,  there  was  remarkably  little  spread  in  the  data  in 
all  these  latitude  bins.  In  the  southern  midlatitudes  (Figure  16),  there  is  less  variability  in 
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Figure  12.  Average  limb  radiance,  and  standard  deviation,  for  SABER  channel  1  for  all  events  in  the 
tirade  bin  30-50°N  on  March  3,  2002;  and  the  same  for  channel  7  nighttime  events.  The  solid  lines  give 
ie  radiance  envelope,  e  g.  the  minimum  and  maximum  radiance  from  any  event  on  that  day. 
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Figure  13.  As  Figure  12,  but  for  the  latitude  bin  10-30°N. 
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Figure  14.  As  Figure  12,  bur  for  the  equatorial  latitude  bin  10°S-10°N. 
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Figure  15.  As  Figure  12,  but  for  the  latitude  bin  10-30°S. 
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C02 15  &  4.3  nm  Mean  Radiance,  30-50S 


Figure  16.  As  Figure  12,  but  for  the  latitude  bin  30-50°S. 


the  mesosphere.  However,  the  4.3  pm  envelope  shows  a  distinct  skewing  toward  large 
radiance  values  at  tangent  heights  above  90  km.  This  is  due  to  one  or  possibly  more 
events  (but  not  many  events)  that  viewed  regions  with  auroral  excitation. 

In  the  high  northern  latitudes  (Figures  17  and  18),  there  is  ample  evidence  of  auroral 
excitation.  Not  only  is  the  envelope  greatly  skewed  to  high  radiance,  but  the  standard 
deviations  fall  outside  the  envelope  on  the  low-radiance  side.  This  is  a  consequence  of 
averaging  over  two  distinct  populations,  which  makes  the  quantities  plotted  less  meaning¬ 
ful  than  they  would  otherwise  be.  The  continued  need  to  look  at  high  latitudes  makes  it 
desirable  to  find  some  way  of  identifying  events  with  auroral  influences  among  the 
SABER  data,  so  that  statistics  could  be  calculated  for  quiescent  events  alone  and  occur¬ 
rence  probabilities  for  auroras  of  different  apparent  strength  could  also  be  determined. 
The  best  way  to  do  that  would  be  with  correlative  data  from  UV  or  visible  sensors,  but 
the  only  instrument  providing  such  data  consistently  is  the  Global  Ultraviolet  Imager 
(GUVI)  on  board  TIMED.  Since  GUVFs  views  are  not  coincident  with  those  of  SABER 
most  of  the  time,  and  since  in  any  case  no  algorithm  has  been  developed  for  accomplish¬ 
ing  this,  the  use  of  SABER  data  alone  is  preferable.  A  very  crude  attempt  to  develop  an 
auroral  discriminant  based  on  channel  7  data  is  described  below. 

It  is  worth  noting,  in  Figures  17  and  18,  the  large  spread  in  4.3  pm  radiance  at  tangent 
heights  of -40-55  km.  This  was  not  seen  in  lower  latitude  regions  on  this  day.  It  is  likely 
due  to  variability  in  the  stratopause  temperature,  but  it  is  interesting  that  it  is  not  apparent 
in  the  15  pm  data. 
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Figure  17.  As  Figure  12,  but  for  the  latitude  bin  50-65°N.  Auroral  effects  skew  the  channel-7  statistics,  as 
discussed  in  the  text. 
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Figure  18.  As  Figure  17,  but  for  the  latitude  bin  65-82°N. 
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Figure  19.  Latitude  distribution  of  channel-7  radiance  at  110  km,  for  all 
nighttime  events  of  March  3,  2002. 


Earlier,  we  mention¬ 
ed  the  need  for  an 
indicator  of  auroral  do¬ 
sing.  Figure  19  shows 
the  distribution  in  lati¬ 
tude  of  channel-7  radi¬ 
ance  at  110  km  for 
March  3.  One  can  see 
greatly  elevated  radi¬ 
ance  levels  for  events 
poleward  of  50  degrees 
in  each  hemisphere.  We 
have  made  use  of  these 
data  to  develop  an  ad 
hoc  auroral  discrimin¬ 
ant.  The  procedure  con¬ 
sists  of  least-squares 
fitting  the  low-latitude 
data  to  a  tenth-order 


polynomial,  determining  the  standard  deviation  of  these  points  with  respect  to  the  fit,  and 
then  asserting  that  any  high-latitude  event  with  radiance  greater  than  five  standard  devia¬ 
tions  beyond  an  extrapolated  baseline  is  an  “auroral  event”.  The  reason  for  using  a 
polynomial  can  be  seen  in  the  shape  of  the  data  near  the  equator,  where  there  is  a  small 


Figure  20.  Channel-7  average  and  standard  deviations,  distinguishing  between  auroral  and  non- 
auroral  events. 
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but  significant  variation  with  latitude.  We  also  note  that  the  extrapolation  to  high  latitudes 
simply  uses  a  constant  value  at  each  end.  And  finally,  deviations  from  fits  from  several 
different  tangent  heights  above  100  km,  not  just  the  one  shown,  were  taken  for  each  event 
in  an  attempt  to  get  a  less  arbitrary  index.  Nevertheless,  a  considerable  degree  of  arbitrar¬ 
iness  does  still  attach  itself  to  the  procedure. 

Figure  20  shows  that,  for  channel  7  at  least,  two  distinct  populations  of  events  can  be 
distinguished  by  our  ad  hoc  approach.  Considering  what  was  done,  that  would  certainly 
be  expected.  In  Figure  21  we  show,  however,  that  other  SABER  data  are  correlated  with 
the  “strength”  of  an  auroral  event,  defined  by  a  number  of  standard  deviations  by  which  a 
given  event’s  channel-7  radiance  exceeds  the  high-latitude  baseline  at  several  tangent 
heights.  This  plot  gives  channel-8  data  for  a  day  in  April  2002  when  a  solar  storm  was  in 
progress  and  when,  therefore,  there  were  many  more  high-latitude  events  with  elevated 
radiance.  We  do  not  plot  the  one-sigma  bars  for  the  “moderate”  and  “strong”  auroral 
events  because  they  are  so  large  that  they  overlap  and  would  obscure  each  other.  The  fact 
that  they  are  so  large  is  suggests  that  the  correlation  between  “strong”  auroral  signatures 
in  channels  7  and  8  is  poor,  and  likewise  between  “moderate”  ones.  This  is  what  would 
be  expected  for  localized  sources,  in  consideration  of  the  fact  that  one  emission  is  optical¬ 
ly  thick  and  the  other  is  not,  and  also  of  the  fact  that  the  index  we  have  defined  has  only 
the  most  tenuous  connection  with  physical  processes  that  produce  strong  and  weak  aur¬ 
ora.  Nonetheless,  one  can  see  from  Figure  21  that  there  is  at  least  some  utility  in  trying  to 
make  the  distinction.  One  point  is  that  the  entire  SABER  data  set  can  be  examined  with 
this  discriminant,  without  resort  to  incomplete  or  hard-to-obtain  ancillary  data. 


SABER  Limb  Radiance,  Channel  8,  50°-65°S 


Figure  21.  OH  data  from  an  active  day  in  April  2002,  distinguishing  quiescent  events  from  those  viewing 
“moderate”  and  “strong”  auroras  according  to  the  index  discussed  in  the  text. 
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We  investigated  one  other  aspect  of  the  March  3  data  that  was  only  briefly  touched 
upon  in  the  earlier  discussion.  The  peaks  and  slope  changes  in  the  15  pm  radiance  pro¬ 
files  that  are  seen  in  many  of  the  equatorial  events  on  that  day  (as  in  Figures  9  and  10) 
suggest  that  temperature  inversions  are  ubiquitous  in  the  upper  mesosphere  in  this  geo¬ 
graphical  region.  Moreover,  the  peaks  typically  appear  at  different  tangent  heights  in  day¬ 
time  and  nighttime.  To  get  a  more  complete  picture,  we  averaged  daytime  and  nighttime 
scans  separately  and  compared  the  averages.  The  results  for  the  equatorial  latitude  bin  are 
given  in  Figure  22,  in  a  fashion  similar  to  that  of  Figures  12-18.  Note  that  the  nighttime 
error  bars  have  small  caps,  whereas  (in  order  to  avoid  cluttering  the  plot)  the  daytime 
ones  do  not.  A  plot  of  the  day/night  difference  as  a  percent  of  the  daily  mean  (including 
night  and  day)  is  also  included.  The  local  time  at  the  ascending  and  descending  nodes 
was  approximately  1700  and  0200,  respectively,  for  all  orbits. 

The  day/night  differences  in  Figure  22  are  striking.  Since  the  15  pm  radiance  is  deter¬ 
mined  primarily  by  the  thermal  structure  of  the  atmosphere,  one  must  conclude  that  the 
wavelike  radiance  pattern,  with  the  day  and  night  profiles  apparently  just  out  of  phase 
with  each  other,  is  a  manifestation  of  the  diurnal  tide.  However,  confirmation  of  that 
conclusion  requires  detailed  modeling,  which  we  have  not  done.  When  reliable  retrieved 
temperature  profiles  become  available,  it  should  be  possible  to  be  more  definite. 


Figure  22.  (a)  SABER  15  pm  limb  radiance  in  the  latitude  bin  spanning  the  equator  for  March  3,  separately 
averaged  over  all  daytime  and  nighttime  events.  Error  bars  indicate  one  standard  deviation  in  each  case,  (b) 
percent  difference  (day  minus  night)  expressed  as  a  percentage  of  the  total  daily  mean  radiance. 
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SABER  Channel  1  Day-Night  Comparison,  30-50N  Day/Night  Difference 


Figure  23.  As  Figure  22,  but  for  the  northern  midlatitude  region. 


SABER  Channel  1  Day-Night  Comparison,  30-50S  Day/Night  Difference 


Figure  24.  As  Figure  22,  but  for  the  southern  midlatitude  region. 
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Figures  23  and  24  present  similar  results  for  northern  and  southern  midlatitude 
regions,  respectively.  For  midlatitudes  one  expects  the  semidiurnal  tide  to  dominate,  and 
indeed  the  day/night  differences  do  not  resemble  the  structures  seen  for  the  equatorial 
region  in  Figure  22.  However,  we  reiterate  that  the  role  that  tidal  processes  may  play  in 
determining  the  15  pm  radiance  has  not  been  investigated  in  detail,  and  comments  relat¬ 
ing  to  it  remain  quite  speculative. 

It  is  interesting  to  note  in  Figures  23  and  24  that  the  nighttime  average  radiance  ex¬ 
ceeds  the  daytime  average  radiance  near  the  mesopause.  All  other  things  being  equal,  one 
expects  non-LTE  effects  to  produce  exactly  the  opposite  result  in  this  region.  In  the  day¬ 
time,  energy  transfer  from  solar-pumped  4.3  and  2.7  pm  states  to  the  bend-stretch  mam- 
fold  elevates  vibrational  temperatures  of  some  of  the  15  pm  states  by  ~4-7  K  for  typical 
midlatitude  conditions,  compared  to  nighttime.  This  produces  an  increase  in  limb  radi¬ 
ance  of  ~5%.  (For  cold  polar  summer  conditions,  much  larger  effects  would  be  expected.) 
Therefore,  the  SABER  results  show  that  all  other  things  are  not  equal,  moreover  that  the 
underlying  conditions  of  the  atmosphere  counteract  and  overwhelm  the  non-LTE  effects. 

The  figures  presented  in  this  section  barely  scratch  the  surface  of  the  information  that 
will  ultimately  reside  in  the  SABER  database.  The  radiance  data  alone  will  provide  a 
very  complete  picture  of  global  radiance  variability  in  much  of  the  infrared  spectral  reg¬ 
ion.  Quantities  that  will  be  retrieved  from  them  (such  as  temperature,  minor  constituent 
densities,  cooling  rates,  etc.)  will  constitute  an  additional  unprecedented  and  very  rich 
resource  for  research  on  the  MLT  region. 
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2.  Atmospheric  Radiance  Code  (ARC) 

During  the  development  of  its  temperature  and  constituent-density  retrieval  algori¬ 
thms,  the  SABER  project  required  a  considerable  amount  of  detailed  non-LTE  modeling. 
This  was  done  to  insure  that  the  operational  code  would  incorporate  the  best  forward 
calculations  that  could  be  performed,  and  to  determine  systematic  and  random  errors  that 
could  result  from  uncertain  input  parameters  (for  example,  rate  constants)  or  imprecise 
numerical  procedures.  Benchmark  calculations  for  many  CO2  states  were  performed 
using  our  non-LTE  model,  Atmospheric  Radiance  Code  (ARC);  some  of  them  were  dis¬ 
cussed  by  Wintersteiner  [2001].  The  need  for  accurate  results  prompted  us  to  write  and 
check  two  new  codes,  which  are  described  very  briefly  below. 

2.1  RADV 

RADV  is  a  new  component  of  ARC.  It  replaces  much  of  the  functionality  of  an  older 
code,  VPMP.  Its  purpose,  like  that  of  RAD  and  RADC  [Wintersteiner  et  al,  1996],  is  to 
calculate  vibrational  temperatures  of  certain  CO2  states.  In  this  case,  the  states  are  the  so- 
called  “Group  1”  {10012-02211-10011}  and  “Group  2”  {11112-03311-11111}  states, 
which  lie  approximately  3600  and  4300  cm"1  above  ground,  respectively.  They  are 
strongly  coupled  to  each  other,  and  also  to  ground  or  lower-lying  C02  states  by  radiative 
transitions  and  by  fast  V-V  exchange  with  N2. 

The  basic  difference  between  RADV  and  VPMP  is  the  inclusion  of  radiative  transfer 
within  the  atmosphere  in  the  new  code,  leading  to  more  accurate  excitation  rates  for  the 
states  in  question. 

In  the  daytime  the  principal  excitation  mechanism  for  the  states  of  Groups  1  and  2  is 
solar  pumping  at  2.7  pm.  Once  excited,  they  emit  photons  at  both  4.3  and  2.7  pm.  This 
contributes  heavily  to  the  total  emission  at  4.3  pm,  and  even,  indirectly,  to  emission  at  15 
pm  (because  the  lower  states  of  the  4.3  pm  transitions  are  in  the  bend-stretch  manifold). 
They  also  affect  the  population  of  vibrationally-excited  N2,  which  feeds  back  and  in  turn 
influences  the  populations  of  other  CO2  states  that  emit  at  4.3  pm.  The  daytime  problem 
was  discussed  in  considerable  detail  by  Sharma  and  Wintersteiner  [1985],  among  others. 
At  night,  there  is  no  strong  excitation  mechanism,  and  these  states  are  often  ignored. 

The  need  to  model  the  emission  at  4.3  pm  very  accurately  caused  us  to  study  the  effect 
of  radiative  transfer  in  2.7  and  4.3  pm  bands  on  the  populations  of  the  Group  1  and  2 
states.  Previously  neglected  because  of  the.  dominance  of  solar  excitation  in  daytime  and 
the  relative  unimportance  of  these  states  at  night,  this  process  was  included  in  RADV 
using  the  basic  ARC  algorithm  for  radiative  transfer  [Wintersteiner  et  al,  1992].  Each 
group  is  connected  to  lower  states  through  two  2.7  pm  bands  and  three  4.3  pm  bands.  The 
results  show  that  the  latter  are  considerably  more  important  than  the  former  for  exciting 
the  high-lying  states. 

In  the  daytime,  the  inclusion  of  this  additional  excitation  mechanism  increases  the  vib¬ 
rational  temperatures  of  the  three  Group- 1  states  by  approximately  2  K  in  the  lower 
mesosphere,  but  only  for  the  major  (626)  isotope  as  shown  in  Figure  25.  This  produces  a 
minor  but  not  negligible  enhancement  in  the  total  4.3  pm  limb  radiance  near  65  km.  For 
the  minor-isotope  and  Group  2  states,  the  changes  amount  to  less  than  half  a  degree  in  the 
mesosphere.  Because  these  states  have  much  lower  absolute  densities,  the  resulting  effect 
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Group  1  Vibrational  Temperature  Comparison: 

Include  &  exclude  radiative  transfer  Effect  of  Including  RT 


Figure  25  (a)  Typical  daytime  vibrational  temperatures  of  the  CO2  major-isotope  (626)  Group  1  states, 
calculated  with  and  without  radiative  excitation  from  within  the  atmosphere,  (b)  Differences  in  vibrational 
temperature  between  the  two  cases. 


on  the  radiance  is  extremely  small.  The  increases  that  are  seen  at  much  higher  altitudes, 
in  the  thermosphere,  also  have  negligible  impact. 

For  the  nighttime  calculation,  in  contrast,  the  inclusion  of  radiative  transfer  is  neces¬ 
sary  for  a  realistic  result.  Whereas  the  Group- 1  major-isotope  states  add  only  a  few  per¬ 
cent  to  the  total  4.3  pm  limb  radiance  at  tangent  heights  in  the  lower  mesosphere,  it  is  still 
worthwhile  to  have  the  capability  to  figure  their  contribution. 

RADV  replaces  VPMP  for  purposes  of  calculating  the  populations  of  Group  1  and  2 
states,  but  the  latter  code  is  still  used  for  the  higher-lying  Group  3  states.  In  daytime, 
these  are  solar-pumped  at  2.0  pm  as  well  as  at  2.7  and  4.3  pm,  and  the  radiative-transfer 
contribution  is  likely  to  make  very  little  difference  Their  significance  at  night  is  com¬ 
pletely  nil. 

2.2  OH_MODEL 

OH  MODEL  is  another  new  component  of  ARC.  It  performs  a  very  simple  calcula¬ 
tion  of  the  populations  of  vibrational  states  of  the  OH  ground  electronic  state,  following 
chemical  production  by  the  so-called  “H  plus  O3”  reaction: 
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i/  +  03 ->0i/(v)  +  O2  v<9 

The  motivation  for  developing  this  code  was  threefold:  to  model  the  SABER  limb 
radiance  in  channels  8  and  9,  which  capture  Av=2  emissions  from  OH;  to  model  the 
contribution  of  Av=l  emission  from  OH  in  SABER  channel  7;  and  to  calculate  the  effect 
that  OH  may  have  on  the  N2  and  CO2  nighttime  vibrational  populations  through  V-V 
transfer,  according  to  a  long-standing  hypothesis  of  Kumer  [ 1978] .  Work  making  use  of 
these  capabilities  is  ongoing  but  incomplete.  We  also  used  OH_MODEL  to  check  results 
from  the  SABER  operational  code,  by  running  cases  with  the  same  input. 

The  new  code  requires  that  one  specify  the  input  model  atmosphere  plus  the  density 
profiles  of  both  atomic  hydrogen  and  ozone.  It  determines  the  rate  of  production  of  the 
nascent  states  (v=6  through  v=9),  and  the  production  and  loss  rates  of  every  level  due  to 
collisional  and  radiative  processes,  and  performs  a  steady-state  calculation  by  equating 
production  and  loss.  The  quenching  scheme  follows  the  prescription  of  Adler-Golden 
[1997]  and  the  Einstein  coefficients  are  from  Goldman  et  al  [1998].  The  correct  values 
for  all  the  parameters  involved  is  still  a  matter  of  some  uncertainty,  and  indeed  has  been 
explored  by  others  at  considerable  length.  However,  we  did  not  do  any  systematic  study 
along  those  lines;  instead  we  relied  on  what  had  already  been  done. 

The  results  that  we  obtained  with  OH_MODEL  are  entirely  as  expected.  The  altitude 
of  the  OH  layer  that  it  predicts  depends  on  the  input  profiles  of  H  and  O3,  as  do  the  num¬ 
ber  densities  that  it  calculates.  The  relative  state-to-state  populations  depend  primarily  on 
the  quenching  scheme  and  the  Einstein  coefficients.  There  was  excellent  agreement  with 
the  SABER  operational  code,  as  anticipated.  The  code  is  now  an  accepted  component  of 
ARC,  and  is  included  as  a  deliverable  to  AFRL,  as  indicated  in  Appendix  A. 
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3.  Summary 

During  the  period  covered  by  this  contract,  we  have  made  considerable  technical 
progress  on  several  fronts  in  fields  that  are  of  interest  to  the  Air  Force.  The  major  ad¬ 
vances  are  detailed  in  the  first  technical  section  of  this  report,  Section  1 . 

We  have  conducted  two  rather  different  investigations  of  atmospheric  radiance  clutter. 
In  the  first,  we  studied  atmospheric  gravity  waves  from  a  physics  standpoint  in  order  to 
develop  a  capability  for  predicting  upwelling  wave-generated  radiance  in  the  CO2  4.3  pm 
band.  There  were  several  motivations  for  doing  this.  One  was  to  replicate  MSX  data 
obtained  in  1 996  with  downlooking  instruments  and  thereby  improve  our  understanding 
of  processes  involved,  including  non-LTE  processes.  Another  was  to  model  the  condi¬ 
tions  and  viewing  geometry  for  which  waves  could  be  expected  to  be  seen,  in  general. 
Although  we  did  not  discuss  it  explicitly  in  our  report,  a  third  motivation  was  to  guide  file 
planning  for  a  proposed  satellite  mission  dedicated  to  studying  atmospheric  gravity 
waves,  Waves  Explorer. 

The  other  investigation  took  the  rather  different  tack  of  examining  new  infrared  radi¬ 
ance  data  obtained  by  the  SABER  instrument  aboard  the  TIMED  satellite.  In  this  case, 
the  idea  is  to  look  at  mean  radiance  and  radiance  variability  over  the  globe  for  different 
geophysical  conditions.  The  ultimate  purpose  would  be  to  catalog  these  results  for  the  full 
range  of  latitudes,  seasons,  local  time,  solar  activity,  etc.,  or  at  least  as  much  of  those 
ranges  as  will  eventually  be  represented  in  the  database.  The  work  reported  in  Section  1 .2 
makes  use  of  only  a  single  day’s  worth  of  data  from  early  in  the  mission,  but  it  shows 
results  that  are  (presumably)  typical  near  equinox.  We  have  given  typical  radiance  pro¬ 
files,  and  daily  or  nightly  means,  and  discussed  their  characteristics.  We  specifically 
pointed  to  effects  caused  by  auroral  activity  and  (we  believe)  tidal  perturbations. 

The  more  important  codes  that  we  have  developed  in  the  course  of  our  investigations 
have  been  made  available  to  researchers  at  the  Air  Force  Research  Laboratory,  along  with 
supporting  files  and  documents.  They  have  been  delivered  in  the  form  of  source  modules, 
as  outlined  in  Appendix  A.  The  two  new  components  of  ARC  are  described  briefly  in 
Section  2.  We  have  also  engaged  in  discussions  about  how  best  to  use  all  the  codes,  and 
what  further  development  of  them  would  be  most  useful  in  the  future. 
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APPENDIX  A:  ARC  AND  SABER  COMPUTER  CODES 

This  Appendix  gives  a  brief  description  of  the  contents  of  the  ZIP  disk  that  has  been 
delivered  to  satisfy  one  of  the  Contract  Data  Requirements,  “computer  software  product 
end  items”.  The  disk  was  written  under  LINUX  but  should  be  readable  by  Windows 
machines  as  well. 

A.1  ARC 

The  current  version  of  Atmospheric  Radiance  Code  (ARC)  can  be  found  in  the  subdir¬ 
ectories  of  the  first-level  directory,  ARC.  As  indicated  in  the  table  below,  there  are  eight 
such  subdirectories.  These  correspond  to  the  seven  primary  codes  that  comprise  ARC, 
plus  the  convolution  program  CONV.  Each  of  the  first  seven  contains  all  the  source  code 
that  is  necessary  for  running  that  program,  as  well  as  a  makefile  that  can  be  used  to  pro¬ 
duce  an  executable.  There  are  between  three  and  seven  source  files  in  these  subdirector¬ 
ies,  each  of  them  containing  one  or  more  concatenated  source  modules.  The  name  of  the 
subdirectory  indicates  what  the  current  version  number  is  for  each  code. 

Each  main  ARC  subdirectory  also  contains 
examples  of  so-called  directives  files.  These  are 
used  to  specify  how  the  respective  programs 
operate:  which  options  to  invoke,  which  input 
files  to  use,  what  geophysical  information  to 
use,  what  values  to  use  for  each  program  para¬ 
meter,  and  what  names  to  give  the  output  files. 
These  are  referred  to  as  “unit-1  files”  because 
they  are  associated  with  logical  unit  1  in  each 
code,  and  they  have  names  with  the  program 
name  and  the  mnemonic  “001”  in  them.  For 
example,  “rad001_nu2_n.dat”  is  a  sample  dir¬ 
ectives  file  for  the  nighttime  CChto)  calcula¬ 
tion  performed  by  RAD. 

The  CONV  subdirectory  contains  a  single 
source  module,  and  a  defaults  file  that  is  open¬ 
ed,  updated,  and  saved  each  time  the  program  is 
run. 

Six  of  the  ARC  codes  were  described  in  detail  by  Wintersteiner  et  al  [1996].  The  other 
two,  RADV  and  OH_MODEL,  are  discussed  in  Section  2  of  this  report. 

A.2  SABER 

The  SABER  directory  contains  two  codes  that  we  have  used  to  process  SABER  Level 
IB  data.  Both  codes  are  works  in  progress,  but  have  been  used  successfully  in  recent 
evaluations  of  SABER  data.  A  complete  description  of  each  code  will  be  prepared  at  such 
time  as  the  code  development  is  complete. 

Briefly,  the  purpose  of  READ_SDAT  is  to  enable  researchers  to  access  and  process 
dozens  or  hundreds  of  LIB  radiance  data  files,  each  corresponding  to  a  single  SABER 
“event”.  The  purpose  is  to  perform  various  tasks  such  as  averaging,  calculating  statistics, 


Table  A-l :  ZIP  Directory  Structure 

/ 

ARC 

RAD_53 

RADCJ7 

RADV_3 

SABS_15 

VPMP_10 

OHMOD_l 

NLTEA_13 

CONV_9 

SABER 

READ_SDAT_6 

RAD_DIST_5 
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sorting  according  to  various  criteria,  and  writing  specific  quantities  or  profiles  to  output 
files  for  further  consideration  (e.g.,  plotting  or  analysis).  The  criteria  currently  in  place 
for  sorting  include  the  latitude  of  the  event,  the  local  time,  the  solar  zenith  angle,  and  the 
tangent  height  (or  range  of  tangent  heights)  of  interest,  and  this  can  be  done  for  one  or 
several  of  the  10  radiance  channels.  There  is  also  a  tentative  indicator  for  the  strength  of 
auroral  activity  seen  by  the  channel-7  detector  that  can  be  used  for  sorting. 

RAD_DIST  is  a  code  that  performs  two  specific  tasks  using,  as  input,  a  particular  type 
of  file  written  by  READ_SDAT.  One  purpose  is  to  sort  events  according  to  the  radiance 
levels  at  particular  tangent  heights,  and  to  create  output  files  giving  the  distribution  of 
these  levels  (e.g.,  the  number  of  events  with  radiance  levels  within  each  bin).  Another  is 
to  generate  the  auroral  indicator  mentioned  above. 

Each  of  the  subdirectories  contains  two  files  with  source  code.  In  each  case,  one  file 
contains  the  main  program  and  the  other  contains  all  the  remaining  source  modules  that 
are  needed.  In  addition,  each  program  has  an  associated  directives  file.  Samples  of  these 
can  be  found  in  each  subdirectory.  These  directives  files  have  the  names  SDAT_NAMES 
and  DIST_NAMES,  and  they  each  supply  specific  information  needed  to  run  their  codes. 

Very  specific  details  about  the  tasks  performed  by  each  of  these  codes,  the  require¬ 
ments  for  performing  them,  and  the  quantities  to  be  read  from  the  directives  files  are 
given  at  the  beginning  of  the  main  source  modules. 


